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I . PURPOSE 


The  purpose  of  this  report  is  to  describe  a calculational  method 
which  can  be  employed  to  provide  a quantitative  measure  of  the  increased 
survivability  of  personnel  in  combat  vehicles  against  the  external  radi- 
ation exposure  threat  posed  by  residual  radiation.  Comparisons  with 
certain  applicable  experiments  are  made  in  order  to  assess  the  accuracy 
of  the  technique. 


ir.  INTRODUCTION 

A.  Background 

Residual  radiation  can  pose  a military  threat.  Under  certain 
meteorological  and/or  weapon  employment  conditions,  radioactive  fallout 
from  a weapon  burst  could  be  of  such  a localized  character  so  as  to 
cause  areas  of  intense  radiation.  The  extent  and  intensities  of  these 
radiation  fields  can  be  mapped  by  troops  utilizing  the  appropriate  equip- 
ment. Alternately,  with  some  information  on  the  conditions  pertaining 
to  the  weapon  burst  and  the  prevailing  winds,  the  fallout  patterns  can 
be  predicted  by  the  battle  field  commander. ^ Troops  passing  through  or 
emplaced  in  these  areas  would  be  subjected  to  a radiation  hazard.  Their 
vulnerability  can  be  reduced  by  a number  of  techniques  which  include 
physically  removing  the  radioactive  debris  by  scraping  the  surrounding 
ground  or  by  utilizing  the  shielding  afforded  by  armored  vehicles.  The 
concern  here  is  the  reduction  of  the  fallout  threat  by  armored  vehicles. 

The  protection  afforded  by  a combat  vehicle  against  the  gamma 
radiation  from  fallout  can  be  quantified  in  terms  of  a gamma  protection 
factor  CGPF)  . The  GPF  for  fallout  radiation  is  defined  as  the  ratio  of 
1^-.  the  dose  at  an  altitude  of  3 feet  in  the  absence  of  the  vehicle  (free 

field  dose)  to  the  dose  at  a position  inside  the  vehicle.  Clearly,  the 
larger  the  overall  GPF  for  a particular  vehicle,  the  more  protection  it 
affords . 

Heretofore,  GPF  values  for  these  geometrically  complicated  military 
vehicles  have  been  measured  experimentally.  One  technique  used  to  simu- 
late the  fallout  threat  against  a variety  of  vehicles  was  a circulating 
point  source  of  radiation. 2 These  experiments  modelled  an  envisioned 
scenario  of  an  infinite-in-extent  flat  air/ground  interface  with  the 
fallout  source  distributed  uniformly  on  the  ground.  Two  convenient 

3 J.C.  Maloney  and  W.J.  Klerm,  "Department  of  Defense  Land  Fallout 

Prediction  System" j May  197Sj  BEL  Report  Do.  2783,  Ballistic  Research 
Laboratory,  APG~EA,  MD.  (AD  #B004148L) 

2 M.A.  Schmoke  and  W.J.  Post,  "Residual  Radiation  Shielding 

Characteristics  of  the  M60A2E2  Tank",  October  2973,  BRL  Report 
No.  1678,  Ballistic  Research  Laboratory,  APG-AA,  MD.  (AD  #§146731) 
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isotopes,  Co  and  Cs,  were  available  to  approximate  the  various 
stages  in  the  radioactive  decay  of  the  fallout  material. 3 oOco,  however, 
was  used  preferentially  since  it  represents  the  average  energy  gamma 
radiation  present  in  fallout  at  early  times. 

Oak  Ridge  National  Laboratory  (ORNL)  developed  a three  dimensional 
Monte  Carlo  radiation  transport  code,  MORSE'^,  which  has  the  capability 
of  handling  complicated  geometries  via  a combinatorial  geometry  (CG) 
package.^  MORSE  is  an  "off-the-shelf"  computer  code  which  has  been 
successfully  employed  in  the  cost-effective  solution  of  a number  of 
diverse  radiation  transport  problems  such  as  reactor  design  and  the 
shielding  of  military  vehicles  against  initial  radiation.®  The  success- 
ful application  of  MORSE  to  the  problem  of  determining  the  amount  of 
shielding  afforded  by  combat  vehicles  against  fallout  radiation  could  be 
a viable  economic  alternative.  It  can  be  estimated  that,  relative  to 
experimental  measurements,  monetary  savings  of  about  70%  could  be  real- 
ized in  the  absence  of  other  considerations.  However,  since  the  math- 
ematical models  of  the  vehicles  used  for  these  calculations  will  have 
already  been  constructed  for  initial  radiation  calculations,  the  effec- 
tive savings  realized  are  even  greater. 

B.  Calculational  Approach 

The  scenario  to  be  modelled  in  the  calculations  is  the  combat 
vehicle  positioned  at  the  interface  of  an  infinite-in-extent  air-over- 
ground environment.  The  source  material  (fallout)  is  assumed  to  be 
spread  uniformly  on  the  smooth  interface.  Figure  1 depicts  the  vicinity 
of  the  combat  vehicle.  In  concert  with  previous  experimental  work,  it 
is  assumed  that  the  gamma  ray  spectrum  emitted  by  the  fallout  can  be 
approximated  by  ^^Co. 

A gamma  ray  emitted  by  the  source  distributed  on  the  ground  can 
interact  with  a crew  member  by  following  one  of  four  paths:  (a)  direct. 


^ B.E.  Rexroad  and  M.A.  Sahmokej  "A  Point  Source  Circulating  System  for 
Simulating  Fallout  Garma  Radiation",  December  1964,  NDL-TM-15, 

Nuclear  Defense  Laboratory,  APG-EA,  MD. 

^ M.B,  Emmett,  "The  MORSE  Monte  Carlo  Radiation  Transport  Code  System", 
ORNL-4972,  Oak  Ridge  National  Laboratory,  Oak  Ridge,  TN,  1975. 

5 

C.E.  Bugart,  "The  Truth  About  Combinatorial  Geometry  Input", 

(1972),  Unpublished  SAI  memo. 

W.A.  Rhoades,  "Development  of  a Code  System  for  Determining  Radiation 
Protection  of  Armored  Vehicles" , 1974,  ORNL-TM-4664,  Oak  Ridge 
National  Laboratory,  Oak  Ridge,  TN. 
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Figure  1.  Vehicle  in  a Fallout  Field. 
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(b)  scattered  by  tank  material,  (c)  sea  tering  in  the  air  or  (d) 
scattering  in  the  ground.  Combinations  of  scatterings  such  as  a ground 
scattering  (d)  followed  by  scattering  with  the  vehicle  material  (b)  be- 
fore interacting  with  a crew  member  are  also  possible,  but  usually  are 
less  likely  for  a given  gamma  ray.  Clearly,  because  of  the  possibility 
of  the  contribution  of  air  or  ground  scattered  gamma  radiation  to  deliv- 
ered dose,  the  infinite  extent  of  these  two  media  should  be  taken  into 
account.  Scattering  from  the  vehicle  suggests  that  its  total  mass  as 
well  as  mass  distribution  will  play  an  important  role  in  the  dose  re- 
ceived by  crew  members  inside  the  vehicle. 

Because  of  the  symmetry  of  the  ground  with  respect  to  a gamma  ray 
emitted  by  a source  particle  on  its  surface,  the  infinite  thickness  of 
the  earth  can  be  accounted  for  exactly  by  replacing  the  earth  with  an 
equivalent  planar  source.  This  is  done  by  noting  that  gamma  rays 
emitted  isotropically  into  the  ground  subsequently  leave  the  ground 
isotropically.  Therefore,  the  effect  of  ground  thickness  can  be  taken 
into  account  by  using  a "modified  source"  that  is  isotropic  and  has  the 
same  energy  spectrum  that  would  be  observed  for  ^*^Co  placed  on  earth  of 
infinite  thickness . 

The  infinite  air  as  well  as  the  geometric  detail  of  the  combat 
vehicle  are  modelled  directly  using  combinatorial  geometry.  The  mean 
free  path  of  a gamma  ray  in  air  is,  typically,  about  300m.  The  air 
surrounding  the  vehicle  is  at  least  3000m  in  extent. 

The  detector  (crew  member's  position)  is  a point  while  the  source 
is  a very  large  plane  surface.  If  the  problem  were  to  be  treated  via 
Monte  Carlo  techniques  in  the  straight  forward  manner,  a gamma -ray  would 
be  started  from  the  air/ground  interface  and  tracked  while  its  effect 
on  the  dose  was  estimated.  While  a gamma  ray  emitted  anywhere  on  this 
large  plane  might  have  an  effect  at  the  detector,  those  emitted  near 
the  detector  would  be  expected  to  contribute  more  dose  to  that  detector 
than  those  emitted  farther  away.  This  latter  fact  can  be  used  to  im- 
prove the  efficiency  of  the  calculation  by  essentially  performing  the 
calculation  backward  in  time.  Time-reversed,  or  adjoint,  particles 
are  started  at  the  detector  point  and  tracked  until  they  cross  the 
air/ ground  interface.  If  an  adjoint  particle  crossing  the  interface 
would  have  been  emitted  from  the  surface  of  the  ground  (due  to  fallout 
radiation)  then  the  adjoint  particle  is  recorded  as  having  had  an 
effect  at  the  detector. 


Approaching  the  problem  with  MORSE  in  the  adjoint  mode  automatically 
takes  advantage  of  the  fact  that  the  nearer  fallout  contributes  most  to 
the  delivered  dose.  In  addition,  less  time  is  wasted  tracking  gamma 
rays  that  would  have  little  or  no  effect  on  the  crew  member.  Several 
problem  dependent  MORSE  subroutines  used  for  these  calculations  are 
listed  in  Appendix  A. 
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Tlie  gamma  ray  data,  material  cross-sections,  and  Auxier-Snyder 
fluence-to-dose  conversion  factors  were  teiken  from  the  DNA  37-21 
library. 7 The  atomic  composition  of  the  air  and  soil  are  listed  in 
Table  I. 


III.  RESULTS 

A.  Ground  Scattered  Spectrum 

Gamma  rays  with  an  isotropic  distribution  were  started  from  a 
point  source  of  ^*^Co  on  the  surface  of  a cube  of  soil  one  kilometre  on 
a side.  All  the  gamma  rays  exiting  from  the  surface  of  the  cube  were 
scored  as  to  energy  and  importance  (statistical  weight)  . These  exiting 
gamma  rays  included  those  which  scattered  from  within  the  cube  of  soil. 

Table  II  contains  the  energy  bin  structure  of  the  gamma  rays  from 
the  DNA  37-21  group  set.  Notice  that  the  two  gamma  rays  of  ^®Co,  1.17 
and  1.33  MeV,  are  both  within  energy  group  11.  Figure  2 presents  the 
results  of  the  above  calculation.  In  Figure  2 it  is  seen  that  50%  of 
the  leakage  gamma  radiation  has  an  energy  equal  to  the  primary  energy 
group.  This  would  be  expected  since  half  of  the  gamma  rays  are  emitted 
away  from  the  soil  cube.  Furthermore,  the  secondary  peak  at  approxi- 
mately 200  keV  is  due  to  backscattered  radiation  (large  angle  scatter- 
ing) which  has  an  increased  probability  of  escaping  the  cube  of  soil. 

B . Free  Field  Calculations 

The  spectrum  presented  in  Figure  2 appears  reasonable;  however, 
another  check  to  determine  whether  it  will  adequately  represent  the 
v-'  leakage  spectrum  for  the  case  of  interest  can  be  performed.  Simulated 

fallout  fields  have  been  studied  by  a number  of  workers.®*®  Reference  9 
contains  semi -empirical  data  for  the  free  field  dose  rate  at  various 
heights  above  a fallout  field  simulated  with  a point  source  of  ^®Co. 

Figure  3 presents  a plot  of  these  data  versus  height.  The  line 
through  these  points  is  intended  solely  as  a guide.  Plotted  also  are 


^ D.E,  Bartinej  J.R.  Knight,  J.V.  Pace  and  E.W.Rouaein,  '-Rvoduotion 

and  Testing  of  the  DR A Few  Groups  Cross  Section  Library",  October  1975, 
ORNL  - TM-4840,  Oak  Ridge  Rational  Laboratory , Oak  Ridge,  TR. 

® Sahumahyk,  et  al.,  "Measurements  of  Gamma  Radiation  ancf  Gamma  Spectra 
versus  Height  Above  a Fallout  Field  Simulated  with  ^^Co",  Rovember 
1965,  RDL-TR-70 , Ruclear  Defense  Laboratory , APG-EA,  MD. 

9 Sahumchyk,  et  al.,  "Scattered  Radiation  (Skyline)  Contribution  to  an 
Open  Basement  Located  in  a Simulated  Fallout  Field",  December  1966, 
RDL-TR~68,  Ruclear  Defense  Laboratory , APG-EA,  MD. 
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Table  I.  Atomic  Composition  of  Air  and  Soil 


Element 


Atomic  Density  (atoms/barn-cm) 


Air 

Oxygen 

1.1229-5* 

Nitrogen 

4.19948-5 

Argon 

2.51482-7 

Ground 

Oxygen 

3.4795Q-2 

Silicon 

1.15967-2 

Aluminum 

4.88019-3 

Hydrogen 

9.75181-3 

* Read  as  1.1229  x 10 
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Table  II . Gamma  lay  Energy  Group  Structure 


Group  Number 
1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


Bead  OB  1.4  X lo'^ 


Upper  Edge 
(eV) 

1.4+7* 

1.0  + 7 

8.0  + 6 

7.0  + 6 

6.0  + 6 

5.0  + 6 

4.0  + 6 

3.0  + 6 

2.5  + 6 

2.0  + 6 

1.5  + 6 
1.0  + 6 

7.0  + 5 

4.5  + 5 

3.0  + 5 

1.5  + 5 

1.0  + 5 

7.0  + 4 

4.5  + 4 

3.0  + 4 
2.Q  + 4 

1.0  + 4 


0.5  6 
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Figure  3.  Dose  Rate  Above  a ^Co  Simulated  Fallout  Field 


the  results  of  15  adjoint  MORSE  calculations  for  the  dose  rate  as  a 
function  of  height  over  an  infinite,  smooth  plane.  Note  that  the  agree- 
ment is  excellent.  The  greatest  difference  (the  point  at  five  feet)  is 
only  5%  which  is  within  the  uncertainties  of  the  MORSE  results  (esti- 
mated at  10%  for  this  case)  . The  quoted  uncertainty  on  the  experimental 
data  is  6%. 

The  agreement  evidenced  above  indicates  that  (1)  the  adjoint  form 
of  MORSE  with  the  DNA  data  set  and  (2)  the  replacement  of  the  ^^Co  plus 
soil  with  an  equivalent  source  are  indeed  valid  for  the  free  field  case. 

C.  In-Vehicle  Calculations  for  Foreign  Medium  Tank:  Vehicle  A 

The  GPF's  have  been  measured  for  a number  of  foreign  and  US 
vehicles.  These  data  were  obtained  under  field  conditions  using  actual 
vehicles  in  a ^®Co-simulated  fallout  field.  The  data  for  a given  vehicle 
were  obtained  as  a function  of  crew  location.  Masonite  manikins,  con- 
structed to  simulate  the  weight  of  the  average  man,  were  placed  at  the 
crew  positions.  Dosimeters  were  placed  at  the  chest,  abdomen,  groin, 
and  back  of  each  of  the  manikins . 

The  GPF's  for  the  four  crew  positions  of  a foreign  medium  tank  (10), 
referred  to  here  as  Vehicle  A,  have  been  measured  using  this  procedure. 

A combinatorial  geometry  description  of  this  vehicle  had  been  construc- 
ted previously  for  initial  radiation  calculations  (11)  . The  description 
of  the  environment  of  the  vehicle  was  modified  to  represent  the  vehicle 
as  situated  on  an  infinite  fallout  field. 

It  should  be  noted  that  the  vehicle  description  was  constructed 
with  careful  attention  to  the  important  aspects  of  the  initial  radia- 
tion scenario  in  which  the  upper  hull  provides  the  most  significant 
shielding.  Thus,  the  lower  portions  of  the  vehicle  were  mathematically 
described  with  some  what  less  detail  than  the  upper  portion  so  as  to 
reduce  the  time  and  cost  of  those  calculations.  The  wheels  and  tracks 
of  Vehicle  A,  for  example,  are  represented  by  two  rectangular  boxes 
located  on  either  side  of  the  vehicle.  The  material  composition  and 
mass  of  these  rectangular  boxes  had  been  adjusted  appropriately  to 
represent  the  average  density  of  the  wheel/track  portion  of  the  vehicle. 
It  was  felt  that  this  degree  of  detail  was  sufficient  for  the  residual 
(fallout)  radiation  shielding  calculations  and  would  not  significantly 
affect  the  final  results.  Therefore,  no  modifications  were  made  to 
the  description  of  the  vehicle  itself. 

The  CG  description  of  the  crew  contained  a torso,  legs  and  thighs 
and  a head.  Except  for  the  head,  which  was  filled  with  air,  the  rest 


For  exact  reference  contact  A.E.  Rainie,  Ballistic  Research 
Laboratory,  APG,  MD. 

For  exact  reference  contact  A.E.  Rainis,  BRL,  APG,  MD. 
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of  the  body  used  the  composite  "man"  material  described  in  Reference  11. 
The  GPF  was  calculated  at  four  locations  in  and  near  crew  members  of 
Vehicle  A;  just  outside  the  center  of  the  chest  (front),  the  center  of 
the  back,  inside  the  air  filled  head,  and  mid  torso.  This  was  done  to 
assess  the  effect  of  the  GPF  on  the  location  of  the  detector  position. 

Figure  4 presents  a comparison  of  the  gamma  energy  spectra  in- 
vehicle  (loader's  position)  and  free  field.  Note  that  the  primary  shape 
differences  are  in  the  low  energy  region  below  200  keV.  This  type  of 
behavior  is  what  might  be  expected  since  the  photoelectric  absorption 
cross  section  for  iron  increases  rapidly  with  decreasing  energy  below 
200  keV.  12 

The  calculated  values  of  the  GPF  for  the  above  locations  and  the 
average  measured  values  are  shown  in  Table  III.  The  uncertainty  limits 
on  the  experimentally  determined  average  GPFs  reflect  both  the  experi- 
mental error  and  the  spread  of  the  values  of  the  GPFs  measured  at  the 
different  locations  on  the  manikin.  The  values  of  the  fractional 
standard  deviation  (fsd)  given  for  the  calculated  values  provide  a cal- 
culational  "figure-of-merit".  That  is,  they  are  a measure  only  for  the 
stochastic  nature  of  the  calculation  and  not  necessarily  an  estimate  of 
the  discrepancy  with  the  "true"  value. 

The  mid-torso  detector  location  is  observed  to  compare  favorably 
with  the  experimental  average  values.  The  largest  observed  difference 
for  these  comparisons  is  less  than  20%,  with  the  calculated  GPFs  always 
within  the  experimental  uncertainty  limits.  Table  III  also  shows  that 
air-detector  calculations  (front,  back  and  head)  generally  lie  close  to 
one  another  but,  with  one  exception,  are  lower  than  the  mid-torso  calcu- 
lations and  the  experimental  average.  However,  with  the  exception  of 
the  driver's  position,  the  experimental  uncertainties  and  the  calculated 
GPFs  with  their  associated  fsd  show  agreement.  Note  that  the  calculated 
GPFs  for  the  mid-torso  location  of  the  gunner's  and  loader's  positions 
have  fsd  values  which  are  twice  that  of  the  calculations  for  the  detec- 
tors in  air.  This  occurs  because  of  the  location  of  the  detector  inside 
the  mid-torso.  From  this  viewpoint,  a detector  position  in  air  is 
desirable.  Of  the  three  in-air  locations  utilized  for  these  calculations, 
the  center  of  the  air  head  is  the  most  easily  identifiable  in  the  com- 
binatorial geometry  description.  Therefore,  both  for  calculational 
convenience  and  ease  of  reproducibility,  the  air  head  location  will  be 
employed  as  the  standard  position  for  future  calculations. 
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D.  In-Vehicle  Calculations  For  Foreign  Medium  Tank:  Vehicle  B 


Calculations  for  another  vehicle  were  performed  to  further  test 
the  correlation  of  an  air  head  detector  location  for  the  calculations. 
The  values  of  the  GPF  for  the  crew  of  this  vehicle  (referred  to  as 
Vehicle  B)  have  also  been  measured  previously . 13  Additionally,  a CG 
description  of  the  vehicle,  previously  used  for  initial  radiation 
calculations,  was  also  available.!^  This  description  was  modified  only 
to  the  extent  described  in  the  previous  section  for  Vehicle  A.  The 
results  of  the  experiments  and  these  calculations  are  presented  in 
Table  IV. 

Note  that  the  average  GPF  for  all  positions  in  the  experiment  and 
the  calculation  agree  with  each  other  despite  the  fact  that  the  experi- 
mental numbers  are  based  on  detectors  placed  about  the  torso  while  the 
the  detector  location  for  the  calculations  is  in  the  head.  Tlie  corre- 
spondence of  these  values  demonstrates  that  the  GPF  for  an  air-head 
detector  is  similiar  to  the  average  GPF  of  detectors  placed  about  the 
body.  This  reinforces  the  use  of  the  air-head  location  for  future 
calculations . 

Examining  the  GPF  values  for  the  individual  crew  positions,  one 
finds  that  the  calculated  values  for  the  driver  position  appear  to  be 
low  when  compared  to  the  experimental  results.  This  may  be  due  to 
limited  detail  of  the  lower  portion  of  the  vehicle.  Because  the 
location  of  the  driver  is  closer  to  the  source  of  radiation  and  the 
"smeared"  wheels  and  track  than  the  other  locations,  the  GPF  for  that 
position  would  be  expected  to  be  the  most  sensitive  to  differences 
between  the  actual  vehicle  and  the  description  employed  for  the 
calculations . 


12  For  exact  reference,  contact  A.E.  Rainis,  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD. 

14  For  exact  reference,  contact  A.E.  Rainis,  Ballistic  Research 
Laboratory , Aberdeen  Proving  Ground,  MD. 
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IV.  CONCLUSIONS 


The  described  calculational  technique  can  be  employed  to  calculate 
values  of  the  GPF  for  the  crew  of  military  vehicles.  Using  the  compari- 
sons between  calculations  and  experiments  on  the  same  vehicles,  a 
conservative  confidence  limit  of  20%  can  be  assigned  to  the  calculated 
values  of  the  GPF.  For  convenience,  the  air-head  is  adopted  as  the 
detector  location  for  future  calculations. 
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APPENDIX  A 

Problem  Dependent  MORSE  Subroutines 

The  computer  code,  MORSE,  used  for  these  calculations  was,  in  the 
main,  the  standard  version  distributed  by  RSIC.**  Several  subroutines 
necessary  for  dose  estimation  are  problem  dependent  and  have  to  be  user 
supplied.  The  subroutines  peculiar  to  this  problem  are  included  in 
Tables  A1-A3, 
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Table  A1 . MORSE  Subroutine  BAfJKR  (Continued) 


Table  A1  MORSE  Subroutine  BANKR  (Continued) 
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Table  A3.  MORSE  Subroutine  SXE 
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APPENDIX  B 


Run  Stream  for  a Typical  Calculation 

A listing  of  the  input  and  controls  cards  for  a typical  problem  is 
presented  for  the  MORSE  Code  as  implemented  on  the  UNIVAC  1108  Computer. 


Appendix  B.  Run  Stream  for  a Typical  Calculation 
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Appendix  B.  Run  Stream  for  a Typical  Calculation  (Continued) 
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Appendix  B.  Run  Stream  for  a Typical  Calculation  (Continued) 
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Appendix  B.  Run  Stream  for  a Typical  Calculation  (Continued) 
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Appendix  B.  Run  Stream  for  a Typical  Calculation  (Continued) 
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